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Ground-state  energies  of  one-  and  two-electron  silicon  dots  in  an  amorphous  silicon  dioxide  matrix 

Davorin  Babic,  Raphael  Tsu,  and  Richard  F.  Greene 
University  of  North  Carolina  at  Charlotte,  Charlotte,  North  Carolina  28223 
(Received  30  October  1991) 

The  one-  and  two-electron  ground-state  energies  of  a  silicon  sphere  embedded  in  an  amorphous  silicon 
dioxide  matrix  are  calculated  as  a  function  of  the  sphere  size.  The  electron-electron  interaction  and  po¬ 
larization  effects  are  treated  by  perturbation;  our  quantum-mechanical  calculation  is  valid  for  small 
spheres  with  radii  between  10  and  40  A.  For  large  spheres,  classical  electrostatics  is  used.  A  universal 
effective  capacitance  is  defined  in  terms  of  the  difference  in  the  ground-state  energies  of  the  (n  -1-1)-  and 
n-electron  cases,  which  agrees  with  the  usual  concept  of  capacitance  in  the  classical  limit. 


I.  INTRODUCTION 

As  pointed  out  by  Hayashi,’  ultrafine  particles  have 
properties  that  suggest  many  scientific  and  technological 
applications.  The  discovery  of  luminescence  in  electro- 
chemically  etched  silicon  is  an  important  scientific  break¬ 
through.^  Recently,  quantum  confinement  of  microcrys¬ 
talline  silicon  (mc-Si)  embedded  in  an  amorphous  silicon 
dioxide  (a-SiOj)  matrix  showed  rich  structure  in  tunnel¬ 
ing  measurements  attributable  to  quantum  confinement 
in  three-dimensions.^  Raman  scattering  from  auantum 
dots  of  Ge  in  a-Si02  shows  the  expected  downshift  from 
the  position  of  bulk  phonons."*  Small  particles  of  GaAs 
(Ref.  5)  and  CdS  (Ref.  6)  reveal  quantum  confinement  in 
optical  absorption.  The  effects  of  charge  accumulation  in 
quantum  confinement  have  been  under  intensive  studies 
involving  quantum  dots,’  as  well  as  quantum-well  struc¬ 
tures.®'*  It  was  pointed  out  that  charge  accumulation 
plays  a  major  role  in  the  energy  states  of  the  confined 
mc-Si.®  We  are  all  familiar  with  the  striking  difference 
between  the  atomic  spectra  of  the  hydrogen  and  helium 
atoms  caused  by  the  presence  of  additional  charge.  The 
physics  is  even  more  complicated  when  the  difference  in 
the  static  dielectric  constant  of  silicon  (12)  and  that  of  the 
a-Si02  (4)  results  in  induced  polarizations.  Specifically,  a 
single  electron  interacts  with  its  induced  polarization  in 
the  oxide.  With  two  electrons  inside  the  silicon,  electrons 
and  induced  polarizations  interact,  resulting  in  a  compli¬ 
cated  picture.  The  present  work  represents  an  attempt  at 
understanding  the  physics  involved.  Eigenstates  of  two 
electrons  in  such  a  system  can  be  analytically  tractable, 
provided  that  spherical  geometry  is  adopted,  and  may 
provide  insight  into  the  behavior  of  actual  quantum  dots. 
The  use  of  the  perturbation  theory  limits  our  results 
essentially  to  the  ground  state. 

The  difference  in  energy  of  a  two-electron  system  from 
a  one-electron  system  may  be  represented  qualitatively  by 
a  capacitance.  Because  of  the  complicated  interactions, 
this  capacitance  is  not  a  constant  representable  only  by 
geometrical  considerations.  In  a  large  quantum  dot  ap¬ 
proaching  the  classical  limit,  the  use  of  a  constant  capaci¬ 
tance'**  should  be  a  good  approximation.  It  will  be  made 
apparent  in  this  work  that  it  is  not  possible  to  calculate 
the  effective  capacitance  of  our  system  in  the  quantum  re¬ 
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gime  with  electrostatics  only. 

The  ground-state  energies  of  a  silicon  sphere  embedded 
in  an  amorphous  silicon  dioxide  matrix  for  one  and  two 
electrons  are  calculated  as  a  function  of  the  size  of  the 
sphere.  Our  results  are  useful  for  an  understanding  of 
transport  measurements,  such  as  resonant  tunneling,  be¬ 
cause  the  voltages  required  to  bring  in  additional  elec¬ 
trons  to  the  silicon  sphere  are  determined  by  the  ground 
states.  Obvio  ^ly,  excited  states  must  be  included  in  a 
proper  treatment  for  an  understanding  of  the  optical 
response. 

II.  MODEL 

The  model  is  taken  as  a  silicon  sphere  of  radius  a,  em¬ 
bedded  in  an  amorphous  silicon  dioxide  matrix.  The 
spherical  geometry  allows  us  to  make  use  of  the  expan¬ 
sion  of  the  Coulomb  potential  in  terms  of  spherical  har¬ 
monics,  resulting  in  a  separation  of  the  radial  and  angular 
components  of  the  wave  functions.  This,  in  turn, 
simplifies  the  evaluation  of  various  integrals  by  the  ortho¬ 
gonality  relationships  among  the  spherical  harmonics. 
The  actual  barrier  height  of  3.2  eV  between  Si  and  a-Si02 
is  taken  to  be  infinite.  This  approximation  avoids  the 
need  to  treat  the  tailing  of  the  wave  functions  into  the  a- 
Si02  region  for  energy  calculations,  but  not  for  tunneling, 
of  course.  The  presence  of  those  tails  necessitates  a  re¬ 
placement  of  the  dielectric  discontinuity  by  a  smooth 
function  in  order  to  avoid  the  singularity  of  the  associat¬ 
ed  polarization  energy."  Certainly,  the  approximation  of 
3.2-eV  barrier  height  by  an  infinite  barrier  height  should 
be  adequate  for  most  lower  states.  Obviously  this  ap¬ 
proximation  does  not  apply  in  the  case  of 
GaAs/Al;tGa,_,As  dots. 

Furthermore,  the  spherical  quantum  box  must  contain 
enough  silicon  atoms  so  that  the  effective-mass  approxi¬ 
mation  and  the  static  values  for  permittivity  are  applic¬ 
able.  These  conditions  are  nearly  satisfied  for  silicon 
spheres  with  radius  in  excess  of  10  A.  Because  of  the 
presence  of  the  longitudinal  and  transverse  masses,  the 
one-electron  states  of  an  isolated  silicon  sphere  require 
the  use  of  a  variational  approach'^  which  would  have 
been  too  complicated  for  us  to  treat  with  inclusion  of  the 
electron-electron  and  polarization  effects.  Therefore,  we 
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have  replaced  in  this  work  the  longitudinal  and  trans¬ 
verse  effective  masses  by  an  isotropic  mass  of  0.26m,. 
The  relative  permittivities  of  Si,  e,,  and  a-SiOj,  Cj*  are 
taken  to  be  12  and  4,  respectively. 

The  Hamiltonian  of  the  system  is,  as  usual,  the  sum  of 
the  kinetic  and  potential  energies.  The  potential  energy 
is  electrostatic  in  nature,  but  its  form  is  substantially 
complicated  by  the  presence  of  two  dielectrics  separated 
by  a  spherical  boundary.  We  treat  this  problem  by  calcu¬ 
lating  electrostatic-potential-energy  terms  via  the  Poisson 
equation  and  then  using  them  in  the  Schrodinger  equa¬ 
tion. 

In  a  typical  resonant-tunneling  experiment,  the  applied 
voltage  measures  the  energy  difference  between  the 
ground  states  of  one-  and  two-electron  systems.  Our  re¬ 
sults  should  therefore  be  useful  in  understanding  trans¬ 
port  measurements.  On  the  other  hand,  optical  measure¬ 
ments  involve  the  occupation  of  excited  states,  which  is 
not  treated  in  this  work.  As  mentioned  earlier,  the  major 
computational  difficulties  are  avoided  with  the  use  of  an 
infinite  barrier.  As  long  as  the  excited  states  for  a  partic¬ 
ular  sphere  size  are  much  below  the  3.2-eV  Si/a-Si02  bar¬ 
rier  height,  excited  states  may  be  treated  within  our 
present  model  with  only  minor  additional  complexity. 

III.  ELECTROS!  ATICS 

Calculation  of  the  electrostatic  energy  terms  follows 
the  outlines  of  the  work  by  Brus  and  by  Bdttcher  on 
dielectrics  in  static  electric  fields.'"  We  cannot  apply  the 
image-potential  construct,  since  we  deal  neither  with  pla¬ 
nar  boundaries  nor  with  perfectly  conducting  spheres. 
Our  chosen  approach  is  to  evaluate  the  electrostatic 
Green’s  function  for  the  problem  at  hand. 

The  electrostatic  potential  far  away  from  the  sphere  is 
set  to  zero,  leading  to  the  following  expressions  for  the 
Green’s  function  inside  the  sphere: 

Gi„{r,r')  =  - - ^ - TT-f  2  ^/'■'^/(cosy)  ,  (3.1) 

4ireoe,lr-r  1  , 

and  outside  the  sphere 

Go^^(t,t')=  J,B,r-^‘+"P,{cosy)  ,  (3.2) 

I 

in  which  r  and  r'  are  the  position  vectors  of  the  field 
point  and  the  charge  point,  respectively,  and  y  is  the  an¬ 
gle  between  these  vectors,  measured  from  an  origin  at  the 
center  of  the  dot. 

The  coefficients  A/  and  B[  are  determined  by  the  elec¬ 
trostatic  boundary  conditions  at  the  Si/a-Si02  interface. 
With  the  use  of  an  infinite  barrier  height,  the  wave  func¬ 
tions  do  not  have  a  tail  extending  outside  of  the  silicon 
sphere,  so  that  the  coefficients  Bi  are  not  needed  in  the 
evaluation  of  matrix  elements.  Consequently,  we  only  ex¬ 
plicitly  exhibit  the  coefficients  Ai  given  by 


and  two-electron  cases.  Even  in  the  one-electron  case,  an 
electrostatic  term  arises  from  the  interaction  of  the  elec¬ 
tron  and  bound  surface  charge  density  at  the  spherical 
boundary  of  two  dielectrics:  an  electron  induces  the 
bound  surface  charge  density  which  generates  electrostat¬ 
ic  potential  at  the  electron’s  position.  Energy  associated 
with  this  term  must  include  a  factor  of  1/2,  since  it  is  a 
self-interaction  term.  Therefore,  this  energy  of  self¬ 
polarization  becomes 


<A.(r)  =  T2 

I 


9^(f,-e2)(/  +  l)r-' 

47reoe'i  [f  2  + ^1  +  ^ ' 


(3.4) 


where  q  denotes  the  electronic  charge. 

In  the  two-electron  case  there  are  four  terms:  self¬ 
polarization  terms  for  each  electron.  Coulomb  interac¬ 
tion,  and  polarization  interaction.  The  polarization  in¬ 
teraction  arises  as  follows:  one  electron  induces  a  bound 
surface  charge  density  which  interacts  with  the  second 
electron.  The  Coulomb  term  has  the  form 


<i>C^T^,T2)  = 


4fl-£ofilr,-r2l 


(3.5) 


and  from  Eqs.  (3.1)  and  (3.3),  the  polarization  term  is 


i2)  = 


2 


^  ^(  e  I  —  £2  )'■  1  f  2^*/  ( cosy ) 
47r£o£,[£2  +  /(£,+£2)]a^'‘''' 


(3.6) 


IV.  QUANTUM-MECHANICAL  CALCULATION 

The  Hamiltonian  for  the  one-electron  case  consists  of 
the  kinetic  energy  for  the  infinite  barrier  potential, 
(F(r)  =  0,  r<ai  oo,  r>a]  and  self-polarization  energy. 
An  exact  analytical  treatment  of  the  Schrodinger  equa¬ 
tion  that  uses  this  Hamiltonian  is  quite  difficult;  we  there¬ 
fore  resort  to  the  perturbation  theory.  We  can  sec  that 
for  sufficiently  small  wells,  the  kinetic  energy  dominates, 
and  the  self-polarization  energy  can  be  taken  into  ac¬ 
count  by  first-order  perturbation. 

The  spherical  Bessel  functions  are  the  solutions  of  the 
zeroth-order  Hamiltonian  that  includes  the  kinetic  ener¬ 
gy  and  infinite  barrier  potential  terms.  Dealing  with  the 
ground  state,  we  are  only  interested  in  the  lowest  eigen¬ 
function'^ 


t(fQ(r)  =  Njo(7rr/a)Yo^{il)  ,  (4.1) 

where  iV  =  a and  /o=0. 050 6606. 

The  self-polarization  energy  term  defined  by 


=  yg(irr/a)rgo(n) 

_  9^(£,-£2)(/  +  l)r^' 

^  ^ _ *  -  _ 

I  87r£o£,[£2-l-/(e,  ■f£2)]a^''''' 


r^drda 


(£l-£2)(/  +  l)r'' 

47r£oe|[£2+/(ei  +  e2nc^^^' 


(3.3) 


contains  a  dimensionless  series. 


(4.2a) 


The  form  of  the  electrostatic  Green’s  functirn  allows 
identification  of  the  potential  energy  terms  in  both  one- 


2 


i+\ 

£,-t-/(£|+£2) 


f'x^'jo{nx}x^dx  , 
•'  0 


(4.2b) 
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which  is  summed  numerically  and  is  equal  to  0.0 IS  16. 
The  final  form  of  the  self-polarization  energy  for  the 
ground  state  is  thus 


Sireofifl 


0.299  . 


(4.3) 


Note  that  the  self-polarization  energy  scales  as  the  in¬ 
verse  of  the  radius  and  is  proportional  to  the  diflFerence  of 
the  dielectric  constants  of  the  two  media. 

The  total  ground-state  energy  of  our  model  system  is 
represented  by  the  following  expression: 

£,  =  144.6/0^-1- 1.44/a  ,  (4.4) 


where  the  energy  is  given  in  eV  and  the  radius  a  in 
angstroms. 

Table  I  shows  ground-state  energies  for  several  sphere 
radii.  Subscripts  1  and  2  of  energy  E  denote  number  of 
electron;  in  the  system.  The  meaning  of  the  superscripts 
is  the  following:  k  for  kinetic  energy,  s  for  self¬ 
polarization,  C  and  p  for  Coulomb  and  polarization  in¬ 
teractions,  respectively,  of  two  electrons,  and  no  super¬ 
script  for  total  ground-state  energy.  The  two-electron  en¬ 
tries  into  the  table  will  be  described  and  evaluated  im¬ 
mediately.  The  last  row  of  the  table,  marked  A,  shows 
the  difference  in  the  ground-state  energies  between  two- 
and  one-electron  systems. 

The  two-electron  Hamiltonian  includes  one-electron 
terms  as  before  and  the  Coulomb  and  polarization  terms 
of  the  two-electron  interaction 

i/  =  ^(V?+V^)-(-K(r,)-(-F(r2)-l-^,(r,) 

2fn 


-  (Wr,  W ,7  Wr,  ¥o< r, ))  , 
:  to  the  following: 

Iq*  J'  f'joiirxOjliirx^) 


(4.6) 


which  is  reducible  to  the  following 


47rcof  ju 


dxj  . 


(4.7) 


The  dimensionless  double  integral  is  calculated  numeri¬ 
cally  and  its  value  is  0.004  585  45.  The  polarization  ma¬ 
trix  element 


9^(c,-e2)(/  +  l)''i''2^/(cosy ) 

47reoe  1  [  ^2  +  /  ( f  1  +  f  2  ^  ^ ' 

(4.8) 


is  substantially  simplified  after  the  application  of  the 
orthogonality  relations  for  the  spherical  harmonics,  be¬ 
cause  all  terms  except  I  =0  vanish.  The  final  value  of  the 
polarization  energy  is 


47reo6,e-)a 


(4.9) 


We  stress  that  self-polarization  energy  contains  contribu¬ 
tions  from  all  Legendre  polynomials,  while  both 
Coulomb  and  polarization  energies  contain  only  the  1=0 
term. 

The  two-electron  ground-state  energy  can  be  written  as 
£2=289.3/0^7.42/0  .  (4.10) 


-l-0,(r2)-f(^c(r„r2)-b<^p(r„r2)  .  (4.5) 

Once  again,  we  choose  to  treat  the  kinetic  energy  for  the 
infinite  barrier  as  the  zeroth  order.  The  other  terms  are 
treated  by  first-order  perturbation  theory.  The  lowest- 
order  spherical  Bessel  function  is  taken  as  the  spatial 
wave  function  for  each  electron  in  the  ground  state.  An- 
tisymmetrization  is  achieved  through  spin  components. 

The  one-electron  terms  of  the  two-electron  ground- 
state  energy  are  the  same  as  in  the  one-electron  case  dis¬ 
cussed  earlier.  We  need  to  evaluate  only  Coulomb  and 
polarization  matrix  elements.  Evaluation  of  the  Coulomb 
matrix  element  proceeds  in  a  manner  similar  to  the  per¬ 
turbation  treatment  of  the  helium  ground  state.  Its 
form  is 


The  kinetic-energy  term  becomes  equal  to  the  other  com¬ 
ponents  of  the  total  energy  at  the  radius  of  39  A,  there¬ 
fore  we  do  not  list  £2  in  Table  I  beyond  40  A.  For  a 
larger  radius,  a  self-consistent  calculation  cannot  be 
avoided.  The  two-electron  ground-state  energy  and  its 
various  components  are  included  in  Table  I  together  with 
one-electron  energy.  The  difference  between  two-electron 
and  one-electron  ground-state  energies  is  also  shown. 

V.  CLASSICAL  CALCULATION 

The  behavior  of  the  system  for  a  very  large  spherical 
well  approaches  its  classical  limit.  The  length  scale  is  set 
up  by  an  electron  coherence  length  that  is  taken  to  be  100 
A.  For  spherical  wells  with  radii  of  100  A  and  larger,  the 


TABLE  I.  Calculated  one-  and  two-electron  ground-state  energies  from  quantum  mechanics. 


a  (A) 

10 

20 

30 

40 

60 

80 

(eV) 

1.446 

0.362 

0.161 

0.091 

0.040 

0.023 

E\  (eV) 

0.144 

0.072 

0.048 

0.036 

0.024 

0.016 

£,  (eV) 

1.59 

0.434 

0.209 

0.121 

0.064 

0.039 

£5  (eV) 

2.893 

0.723 

0.321 

0.182 

Ei  (eV) 

0.288 

0.144 

0.096 

0.072 

E^  (eV) 

0.214 

0.107 

0.071 

0.054 

£f  (eV) 

0.24 

0.12 

0.08 

0.06 

£2  (eV) 

3.635 

1.094 

0.568 

0.368 

A  (eV) 

2.05 

0.66 

0.36 

0.25 
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electron  kinetic  energy  is  estimated  via  the  uncertainty 
principle.  The  estimate  is  of  the  order  of  1  meV,  which  is 
negligible  compared  with  the  electrostatic  contribution. 
Evaluation  of  the  electrostatic  energy  requires  that  one 
specifies  electronic  positions.  A  single  electron  in  the  sys¬ 
tem  has  its  polarization  energy  minimum  at  the  center  of 
the  well.  The  easiest  way  to  calculate  the  energy  in  this 
case  is  via  energy  density.  The  expression  that  needs  to 
be  evaluated  is 


E 


_  1 

2 

^2 

fl 

1  ^1-^2 
2 


l6iT^€Qr* 


dr 


(5.1) 


leading  to  the  following: 


A7r€Q£^a 


(5.2) 


In  order  to  find  the  positions  of  the  two  electrons,  the 
classical  ground-state  energy  in  the  two-electron  case  is 
determined  by  the  minimum  of  the  competing  repulsive 
components.  Coulomb  and  polarization  interactions 
push  the  electrons  apart,  towards  the  well’s  boundary, 
while  the  self-polarization  terms  push  the  electrons  away 
from  the  boundary,  closer  to  each  other.  The  electrons 
are  located  at  the  same  distance  from  the  spherical-well’s 
center  opposite  to  each  other.  The  distance  from  the 
center,  denoted  by  b,  is  obtained  by  minimization  of  the 
total  electrostatic  energy 


4rreoe, 


1 

—  +  \ 
lb  ^ 


+  2 


-62)b“(/-H )(-!)' 
/  a^'  +  '[e2  +  /(e|  +  e2)] 

(gi-f2)b^'(/  +  l) 

|2'  +  '[e2+/(e,+e2)] 


(5.3) 


The  above  expression  can  be  written  as  a  function  of  the 
ratio  between  an  electron’s  distance  from  the  center  and 
the  well’s  radius,  x=b /a  : 


FIG.  1.  Ground-state  energy  of  one-  and  two-electron  sys¬ 
tems  with  respect  to  the  spherical-well  radius.  The  dashed  por¬ 
tion  of  El  represents  an  extrapolation  to  agree  with  the  classical 
calculation  for  large  a. 


between  the  electrostatic  energy  computed  classically  and 
quantum  mechanically.  The  electrostatic  energy  comput¬ 
ed  classically  is  a  mathematically  well-defined  quantity 
for  radii  less  than  100  A,  although  it  loses  physical  mean¬ 
ing  in  this  range.  We  plot  the  ground-state  energy  of 
one-  and  two-electron  systems  with  respect  to  the 
sj)herical-well  radius  in  Fig.  1.  Since  we  restrict  a  to  39 
A  or  less,  for  a  greater  than  30  A  is  obtained  by  a 
smooth-curve  fit  through  E2  at  39  and  120  A  obtained 
classically. 
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The  minimum  with  respect  to  x  is  found  numerically.  Its 
value  is  5.0284  at  x  =0. 594. 

Table  II  shows  values  of  the  one-  and  two-electron 
ground-state  energies  and  their  difl!erence  computed  clas¬ 
sically  for  several  radii  of  the  spherical  well.  The  mean¬ 
ing  of  the  symbols  is  the  same  as  in  Table  I.  We  include 
the  radii  smaller  than  100  A  to  emphasize  the  difference 


VI.  DISCUSSION  AND  CONCLUSIONS 

In  this  work  an  attempt  to  model  three-dimensional 
quantum  confinement  in  Si  nanometer-sized  crystals  sur¬ 
rounded  by  an  fl-Si02  matrix  is  presented.  We  are  in¬ 
terested  in  the  applied  voltage  required  for  tunneling  of 
initially  one  and  then  two  electrons  through  the  quantum 
well.  A  single  electron  tunnels  when  the  applied  voltage 
aligns  the  Fermi  level  of  the  outside  (degenerate)  electron 
reservoir  with  the  one-electron  ground-state  energy  of  the 
well.  Additional  voltage  is  necessary  to  have  two  elec¬ 
trons  tunnel  through  the  well.  This  voltage  is  equal  to 


TABLE  II.  Classically  calculated  one-  and  two-electron  electrostatic  energies. 


a  (A) 

10 

20 

30 

40 

60 

80 

100 

120 

Ef  (eV) 

0.12 

0.06 

0.04 

0.03 

0.02 

0.015 

0.012 

0.01 

£f  (eV) 

0.60 

0.30 

0.20 

0.15 

0.10 

0.075 

0.06 

0.05 

A  (eV) 

0.48 

0.24 

0.16 

0.12 

0.08 

0.06 

0.048 

0.04 
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the  difference  of  the  ground-state  energies  between  two- 
and  one-electron  systems. 

Our  model  is  a  compromise  between  computational 
feasibility  and  the  need  t'^  include  all  physically  relevant 
terms  in  the  energy.  We  evaluate  the  energy  in  two  com¬ 
putationally  convenient  limits:  for  very  small  wells  in 
which  the  kinetic  energy  dominates,  so  that  polarization 
energy  may  be  treated  as  a  perturbation;  and  for  wells 
large  enough  so  that  the  situation  can  be  considered  clas¬ 
sically.  Our  intuitio'’  led  us  to  view  the  problem  of  two 
electrons  in  a  small  well  as  a  cross  between  the  helium- 
atom  and  particle-in-a-box  problems.  The  helium-atom 
model  supports  the  suitability  of  the  perturbation  theory 
to  deal  with  electron  interactions  in  the  ground  state. 
The  particle-in-a-box  approach  permits  the  use  of  the 
infinite-barrier  potential,  thus  allowing  vanishing  bound¬ 
ary  conditions  for  the  electron  wave  function  at  the 
dielectric  boundary  to  be  :  ;t  up. 

In  the  classical  limit,  in  order  to  apply  electrostatics, 
we  need  to  specify  the  positions  of  the  two  electrons; 
therefore,  a  minimization  of  the  total  electrostatic  energy 
is  necessary.  The  two  electrons  are  located  at 
(r  —  b,d  =  0)  and  {r  =  b,d  =  ir).  Consequently,  the 
quantum-mechanical  wave  functions  for  the  two  elec¬ 
trons  cannot  be  represented  by  the  zeroth-order  spherical 
Bessel  function,  even  for  small  well  sizes.  This  fact  clear¬ 
ly  points  out  the  limitation  in  the  app'icability  of  the 
present  perturbation  calculation. 

It  is  important  to  note  that  all  the  electrostatic  terms 
present  have  been  fully  taken  into  account.  Our  results 
clearly  indicate  the  importance  of  the  self-polarization 
and  polarization  interaction  terms,  which  have  compara¬ 
ble  magnitudes  but  apparently  are  somewhat  larger  than 
the  Coulomb  interaction.  Therefore,  the  use  of  a 
voltage-independent  capacitance^  '®  to  deal  with  the 
effects  of  charge  accumulation  does  not  offer  a  clear  phys¬ 
ical  insight.  The  most  physically  correct  definition  of  the 
capacitance  should  correspond  to  what  may  be  measured. 
Consequently,  we  define  an  effective  capacitance  C^g  in 
terms  of  an  increment  of  the  total  stored  energy  in  the 
system.  For  the  particular  case  of  one-  and  two-electron 
systems,  C^g  may  be  defined  by 

=  (6.1) 

This  universal  definition  of  Cj(f  gives  values  for  smaller 
spheres  far  below  what  can  be  calculated  from  classical 


electrostatics.  Since  we  assume  at  120  A  is  adequate¬ 
ly  represented  by  its  classical  value,  at  60  and  30  A, 
is  only  one-half  and  one-third  as  large  as  the  classical 
values,  respectively.  It  is  important  to  note  that  the  usual 
concept  of  a  capacitance  is  meaningful  in  the  classical 
limit. 

Our  calculations  on  a  suitably  defined  model  system 
cover  the  ground-state  energies  of  relatively  small  and 
relatively  large  quantum  wells.  We  believe  that  ground- 
state  energies  of  the  intermediate  quantum  wells  can  be 
obtained  to  a  fair  degree  of  accuracy  by  a  smooth  inter¬ 
polation  scheme.  We  also  believe  that  we  have  set  the 
correct  voltage  scale  as  high  as  1  V  for  the  problem  we 
dealt  with,  which  is  much  larger  than  the  corresponding 
scale  for  GaAs. 

Extension  of  our  results  to  a  larger  number  of  electrons 
and  excited  states  faces  severe  calculational  impediments. 
The  infinite-barrier  approximation  is  less  valid  in  the  case 
of  a  larger  number  of  electrons  than  it  is  for  the  excited 
states  of  one-  and  two-electron  systems.  The  ut.lL-  ion 
of  a  finite  barrier  immediately  introduces  a  gamut  of 
difficulties:  electron  wave  functions  extend  a  tail  across 
the  dielectric  boundary,  which  prohibits  the  use  of  the 
dielectric-constant  discontinuity  at  the  boundary,  leading 
to  a  complicated  electrostatic  problem.  The  difficulties 
are  compounded  by  the  need  to  mix  in  higher  spherical 
Bessel  functions  and  spherical  harmonics  to  describe 
correctly  the  situation  in  the  well.  Everything  points  to  a 
difficult  problem  that  cannot  be  dealt  with  analytically 
but  only  numerically.  Our  work  has  elucidated  the  major 
physics  of  the  problem  without  sacrificing  too  much 
rigor,  while  maintaining  computational  simplicity  and 
feasibility.  As  pointed  out  previously,  the  problem  in¬ 
volving  low-lying  excited  states  should  be  tractable  for 
one-  and  two-electron  systems,  with  use  of  the  infinite 
barrier.  Therefore,  optical  transitions  between  the 
ground  state  and  the  first  excited  state  should  not  be  toe 
much  more  complicated,  when  treated  within  the  frame¬ 
work  of  our  model. 
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